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Stability  and  hydrogen  storage  properties  of  various  metal-decorated  benzene  complexes  have  been 
studied  using  a  first-principles  method.  According  to  our  results,  most  metals  from  the  first  to  the  third 
period  of  the  periodic  table  can  be  strongly  adsorbed  on  the  benzene  surfaces,  except  Be,  Na,  Mg,  K,  and 
Zn.  Among  the  metals  we  studied  Ca  is  the  most  promising  adsorbate  for  hydrogen  storage.  Two  Ca 
atoms  prefer  to  be  isolated  on  the  benzene  surfaces,  one  on  each  side,  and  each  adsorbed  Ca  can  adsorb 
up  to  four  H2  molecules,  thus  yielding  a  H2  uptake  of  9.2  wt%.  The  calculated  binding  energy  is 
0.45  eV/H2,  suitable  for  reversible  hydrogen  storage.  The  adsorbed  Li  and  Sc  can  also  adsorb  a  large 
number  of  H2  molecules  with  suitable  H2  binding  energies,  but  they  are  likely  to  dimerize  on  the  benzene 
surfaces.  Other  metal  adsorbates  we  studied  are  not  suitable  for  hydrogen  storage,  because  they  will 
suffer  from  the  large  H2  binding  energy  or  the  low  hydrogen  storage  capacity. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  success  of  hydrogen  economy  critically  depends  on  the 
discovery  of  hydrogen  storage  materials  meeting  the  requirement  of 
high  gravimetric  and  volumetric  density,  fast  kinetics,  and  favorable 
thermodynamics  [1-4].  In  recent  years,  light  metal-decorated 
carbon-based  materials  [5,6]  such  as  C6o  fullerene  [7-12],  graphene 
[13—15],  small  organic  molecules  [16-18],  and  metal-organic 
frameworks  (MOFs)  [19-24]  have  attracted  much  attention  for 
potential  applications  in  hydrogen  storage.  In  these  systems,  the 
metal  denotes  electrons  to  the  host  materials  and  occurs  in 
a  cationic  form,  and  then  it  can  adsorb  a  large  number  of  H2 
molecules.  Studies  on  these  systems  do  provide  promising  results, 
and  at  the  same  time,  they  encounter  various  problems.  For 
example,  though  light  metal-decorated  C6o  systems  can  adsorb 
a  large  number  of  H2  molecules  [7,9,12],  they  may  suffer  from  the 
small  hydrogen  binding  energy  [7]  or  the  strong  clustering  ability  of 
the  metal  atoms  [10,11]. 

We  also  find  that  when  considering  the  hydrogen  storage 
properties  of  materials  such  as  C6o,  MOFs  and  so  on,  the  six-carbon 
ring  often  plays  important  roles.  Metals  are  usually  adsorbed  above 
the  six-carbon  ring  and  then  adsorb  H2  molecules,  such  as  alkali 
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metals  in  C6o  [25],  Ti  atoms  in  C6o  [10],  Li  atoms  in  MOF-5  [19]  and 
so  on.  From  this  point,  we  think  that  it  will  be  interesting  to 
investigate  the  hydrogen  storage  capability  and  capacity  of  one 
simple  material  containing  the  six-carbon  ring,  and  to  do  this 
benzene  is  the  simplest  selection.  Besides,  we  note  that  3d-tran- 
sition  metal-benzene  complexes  have  been  successfully  synthe¬ 
sized  experimentally  [26].  This  makes  it  practically  interesting  to 
study  the  hydrogen  storage  properties  of  metal-decorated  benzene 
complexes. 

In  this  paper,  we  present  theoretical  studies  on  the  stability  and 
hydrogen  storage  properties  of  various  metal-decorated  benzene 
complexes.  All  the  eighteen  metals  from  the  first  to  the  third  period 
of  the  periodic  table  are  studied  as  adsorbates.  We  find  that  most 
metals  are  strongly  adsorbed  on  the  benzene  surfaces,  except  Be, 
Na,  Mg,  K,  and  Zn.  However,  the  adsorbed  Al,  Ga  and  transition 
metals  from  Ti  to  Cu  are  not  suitable  for  hydrogen  storage,  because 
they  cannot  meet  the  H2  binding  energy  or  the  H2  gravimetric 
density  requirements.  The  adsorbed  Li  and  Sc  can  adsorb  a  large 
number  of  H2  molecules  per  metal  with  a  suitable  H2  binding 
energy,  but  they  will  dimerize  on  the  benzene  surfaces.  Two  Ca 
atoms  prefer  to  be  isolated  on  the  benzene  surfaces,  one  on  each 
side,  and  each  adsorbed  Ca  can  then  adsorb  four  H2  molecules.  The 
binding  energy  is  always  within  a  few  tenths  of  an  eV  per  H2 
molecule,  suitable  for  reversible  hydrogen  storage.  Further  molec¬ 
ular  dynamics  (MD)  simulations  indicate  that  the  Ca-decorated 
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benzene  complexes  are  stable  at  300  °C.  Based  on  these  discussions, 
we  suggest  that  Ca-decorated  benzene  complexes  are  the  most 
promising  media  for  hydrogen  storage  in  this  paper. 

2.  Model  and  method 

The  whole  work  is  based  on  first-principles  spin-polarized 
calculations.  We  use  the  projector  augmented  wave  (PAW)  [27] 
method  as  implemented  in  the  VASP  codes  [28],  where  the  elec¬ 
tronic  exchange  and  correlation  are  approximated  using  the  local 
density  approximation  (LDA)  with  Ceperly  and  Alder  (CA)  func¬ 
tional  [29].  Single  molecular  complexes  are  treated  in 
a  20  x  20  x  20  A  supercell  with  T  k-point  and  an  energy  cutoff  of 
400  eV.  By  increasing  the  supercell  size  and  the  energy  cutoff  we 
find  that  the  energy  difference  is  less  than  0.02  eV,  indicating  good 
convergence.  During  relaxations,  all  the  ions  are  allowed  to  relax 
until  the  force  acting  on  each  ion  is  less  than  0.01  eVA-1.  The 
calculated  C-C  distance  in  the  benzene  ring  is  1.39  A,  in  good 
agreement  with  experimental  reports  [30].  Furthermore,  we 
compare  the  LDA  results  with  those  of  the  generalized  gradient 
approximation  (GGA)  with  Perdew-Wang  91  (PW91)  [31]  and 
Perdew-Burke-Ernzerhof  (PBE)  [32]  functionals.  We  find  that  LDA 
binding  energies  are  slightly  larger  than  those  of  GGA  results.  For 
example,  in  Ca-decorated  benzene  case  the  H2  binding  energy 
calculated  within  LDA  is  about  0.2  eV  larger  than  that  calculated 
within  GGA.  However,  GGA  calculations  using  PW91  and  PBE  give 
almost  the  same  binding  energies. 

We  examine  various  possible  configurations  of  C6H6Mn 
complexes  and  their  corresponding  H2  binding  properties.  Here,  M 
stands  for  any  one  of  the  eighteen  metals  from  the  first  to  the  third 
period  of  the  periodic  table.  We  only  consider  complexes  contain¬ 
ing  one  (n  =  1)  or  two  (n  =  2)  M  atoms  because  more  metal  atoms 
may  lower  the  hydrogen  storage  capacity  and  thus  should  be 
avoided  in  the  syntheses.  Metal  binding  energy  is  calculated  as 
Eb(M)  =  [nE(M)  +  E(C6H6)-E(C6H6Mn)]/n,  where  E(M),  E(C6H6) 
and  E(C6H6Mn)  are  the  energies  of  the  metal  atom,  the  benzene 
molecule,  and  the  CeH6Mn  complex,  respectively.  The  H2 
binding  energy  is  calculated  as  Fb(H2)  =  [F(C6H6Mn)  +  m£(H2)  - 
E(C6H6Mn  +  mH2)]/m.  In  this  formula,  E(H2)  and  E(C6H6Mn  +  mH2) 
are  the  energies  of  the  H2  molecule  and  the  C6H6Mn  +  mH2 
configuration,  respectively,  and  m  is  the  number  of  the  adsorbed  H2 
molecules. 

3.  Results  and  discussions 

3.1  Hydrogen  storage  in  CgHgM  complexes 

We  start  from  n  =  1  case.  To  study  the  hydrogen  storage  prop¬ 
erties  of  various  CeHeM  complexes,  several  points  should  be  noted. 
First,  different  metals  may  have  different  adsorption  sites  on  the 
benzene  surfaces,  and  some  metals  may  not  be  adsorbed  on  the 
surfaces  completely.  Second,  even  if  some  metals  can  be  adsorbed 
on  the  surfaces,  it  does  not  mean  they  can  further  adsorb  H2 
molecules.  Third,  even  if  some  metals  can  satisfy  the  first  and  the 
second  requirements,  the  metal  adsorption  sites  may  be  different 
before  and  after  H2  adsorption.  To  make  these  questions  clear,  we 
first  place  various  metals  to  different  sites  of  the  benzene  surfaces 
and  then  completely  relax  the  structures.  The  initial  sites  include 
the  H-site  (above  the  center  of  the  carbon  ring),  the  A-site  (above 
one  carbon  atom),  the  S-site  (above  the  C-C  single  bond),  and  the 
D-site  (above  the  C=C  double  bond).  Through  the  relaxed  struc¬ 
tures  as  well  as  comparisons  of  metal  binding  energies  at  different 
adsorption  sites,  we  can  judge  which  metals  can  be  adsorbed  on  the 
benzene  surfaces  and  where  they  will  prefer  to  be  adsorbed.  To  the 
stable  C6H6M  complexes,  H2  molecules  are  added  one  by  one 


around  the  metals  until  they  begin  to  escape.  Through  this  process 
we  can  obtain  the  maximum  number  of  adsorbed  H2  molecules 
Nmax( H2)  on  various  complexes  and  the  atomic  geometries.  And 
finally,  we  can  calculate  the  average  H2  binding  energy  and  the  H2 
gravimetric  density  for  every  case.  Metal  binding  properties  on 
benzene  as  well  as  their  H2  binding  properties  for  n  =  1  case  are 
listed  in  Table  1,  where  only  the  configurations  that  can  stably 
adsorb  H2  molecules  have  been  shown. 

When  a  metal  atom  is  moved  toward  benzene,  the  large  elec¬ 
tronegativity  of  benzene  facilitates  the  transfer  of  electrons  from 
the  metal  atom  to  benzene,  resulting  in  an  electrostatic  interaction 
between  them  and  leaving  the  metal  atom  in  a  cationic  form.  It 
should  be  noted  that  the  metal-benzene  interaction  may  be  related 
not  only  to  the  cationic  charge  state  but  also  to  the  metal  atom  size, 
thus  making  the  interaction  or  the  binding  energy  complicated.  The 
main  group  metal  atoms  such  as  Li  and  Ca  trap  H2  molecules  mainly 
by  an  electrostatic  polarization  mechanism,  while  the  transition 
metal  atoms  do  that  mainly  by  the  Kubas  interaction  [33].  The  H2 
molecules  transfer  their  a  electrons  to  the  unoccupied  transition 
metal  d-orbitals,  and  then  the  transition  metal  back-denotes  its 
occupied  d  electrons  to  hydrogen  a*  antibonding  orbitals.  There¬ 
fore,  H2  binding  energy  and  the  maximum  number  of  adsorbed  H2 
molecules  on  the  adsorbate  should  also  be  related  to  the  properties 
of  the  metal  d-orbitals,  such  as  the  d-energy  and  so  on.  From 
Table  1  we  can  see  that  the  binding  energies  of  Na  and  K  are  always 
negative,  while  those  of  Mg  and  Zn  are  always  less  than  0.2  eV/M, 
no  matter  where  they  are  adsorbed.  The  negative  or  small  binding 
energies  indicate  that  these  metals  cannot  form  complexes  with 
benzene  molecules  stably,  so  they  are  not  suitable  for  hydrogen 
storage.  Group  IIIA  elements  Al  and  Ga  are  likely  to  occupy  the 
A-sites,  while  other  metals  are  likely  to  be  adsorbed  on  the  H-sites. 
It  is  obvious  that  all  the  metals  except  Zn  are  positively  charged 
when  forming  complexes  with  benzene  molecules.  Besides,  charge 
on  the  adsorbed  transition  metals  from  Sc  to  Zn  decreases 
from  +1.06e  to  -0.01  e,  which  is  in  agreement  with  the  increased 
electronegativity  from  Sc  to  Zn. 

We  find  that  though  group  IIIA  metals  Al  and  Ga  have  consid¬ 
erable  metal  binding  energies  on  the  A-sites,  they  are  not  able  to 
adsorb  any  H2  molecules.  Besides,  when  trying  to  place  H2  mole¬ 
cules  around  the  adsorbed  Fe  atom,  they  will  either  escape  or  be 
dissociated.  Therefore,  Fe  is  also  not  a  suitable  adsorbate  for 
hydrogen  storage.  The  adsorbed  Li,  Ca,  Sc  and  Ti  can  adsorb  four  H2 


Table  1 

Metal  and  H2  binding  properties  in  n  =  1  (C6H6M)  case. 


Metal 

M-site  a 

Eb(  M) 
(eV) 

M-chargeb 

Nmaxi  H2) 

Eb(  H2) 
(eV) 

H2  weight 
(wt%)c 

Li 

H 

0.54 

+0.77e 

4 

0.22 

8.6 

Be 

H 

0.74 

+1.24e 

- 

- 

- 

Na 

- 

-0.02 

- 

- 

- 

- 

Mg 

H 

0.10 

+0.05e 

- 

- 

- 

Al 

A 

1.15 

+0.80e 

- 

- 

- 

K 

- 

-0.10 

- 

- 

Ca 

H 

0.47 

+1.02e 

4 

0.56 

6.4 

Sc 

H 

2.28 

+1.06e 

4 

0.69 

6.1 

Ti 

H 

2.94 

+1.22e 

4 

1.05 

6.0 

V 

H 

2.78 

+0.96e 

3 

1.21 

4.4 

Cr 

H 

1.12 

+0.81e 

3 

1.46 

4.4 

Mn 

H 

1.11 

+0.67e 

3 

1.08 

4.3 

Fe 

H 

3.01 

+0.55e 

- 

- 

- 

Co 

H 

3.12 

+0.41  e 

2 

1.26 

2.8 

Ni 

D 

3.10 

+0.3  le 

2 

0.92 

2.8 

Cu 

S 

0.77 

+0.27e 

2 

0.60 

2.8 

Zn 

H 

0.17 

-O.Ole 

- 

- 

- 

Ga 

A 

0.82 

+0.02e 

- 

- 

- 

a  The  most  stable  site  of  various  metal  atoms. 
b  Bader  charge  on  various  metal  atoms. 
c  H2  gravimetric  density. 
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molecules  per  metal  at  most,  forming  Ce^M  +  4H2  configurations. 
In  Li,  Ca  and  Sc  cases,  the  H2  gravimetric  densities  decrease  grad¬ 
ually  (see  Table  1 ),  but  they  are  always  higher  than  the  gravimetric 
density  (6  wt%)  of  the  U.  S.  DOE  target  by  2010  [34];  the  H2  binding 
energies  increase  gradually,  but  they  are  always  within  a  few  tenths 
of  an  eV,  suitable  for  reversible  hydrogen  storage.  This  is  not  the 
case  for  Ti  adsorbate,  which  can  also  adsorb  four  H2  molecules 
per  metal,  but  the  binding  energy  is  1.01  eV/H2,  too  large  to  be 
reversible.  We  note  that  the  binding  energies  in  Ti-decorated  C6o, 
graphene,  and  CNT  cases  are  0.55  [11],  0.35  [35],  and  0.56  [35] 
eV/H2,  respectively,  all  smaller  than  the  value  in  this  work. 
However,  the  value  in  an  isolated  Ti+  ion  case  is  1.07  eV/H2  [36], 
very  close  to  that  in  this  work.  This  is  in  good  agreement  with  bader 
analysis  that  the  adsorbed  Ti  carries  a  positive  charge  of  1.22e  (see 
Table  1 ).  CeHeLi  +  4H2  and  CeHeCa  +  4H2  are  geometrically  similar, 
so  we  only  show  the  structure  of  the  later  one  in  Fig.  1(a)  as  an 
example.  We  can  see  from  the  figure  that  after  H2  adsorption,  the 
adsorbate  still  occupies  the  H-site.  The  calculated  C-Ca  distance  is 
2.15  A,  much  smaller  than  sum  of  the  atomic  radii  of  C  and  Ca 
(3.14  A),  indicating  the  formation  of  C-Ca  chemical  bonds.  Besides, 
four  H2  molecules  distribute  symmetrically  around  the  Ca  atom, 
and  the  calculated  H-H  distance  in  the  adsorbed  H2  molecules  is 
0.83  A.  This  is  slightly  larger  than  the  H-H  distance  (0.75  A)  in  an 
isolated  H2  molecule,  indicating  a  physisorption  in  a  molecular 
manner.  The  atomic  geometry  of  CeHeTi  +  4H2  is  shown  in  Fig.  1(b), 
and  that  of  CeH6Sc  +  4H2  is  also  similar.  The  H-H  distance  is  about 
0.85  A  in  Sc  case  and  about  0.95  A  in  Ti  case.  In  order  to  judge 
whether  the  adsorbed  H2  molecules  have  dissociated,  we  show  the 
charge  distributions  of  C6H6Ti  +  4H2  in  Fig.  1(c).  From  the  figure, 
obvious  charge  distributions  can  be  seen  between  the  adsorbed 
hydrogen  atoms  composed  of  the  H2  molecules,  indicating  the 
hydrogen  storage  is  in  a  molecular  manner. 

The  adsorbed  V,  Cr  and  Mn  can  adsorb  three  H2  molecules  per 
metal  as  listed  in  Table  1,  forming  CeHeM  +  3H2  configurations.  The 
three  cases  have  similar  atomic  geometries,  so  we  only  show  that  of 
CeH6Mn  +  3H2  in  Fig.  2(a)  as  an  example.  From  the  figure  we  can 
see  that  the  three  adsorbed  H2  molecules  distribute  symmetrically 
around  the  adsorbate.  The  calculated  H-H  distances  in  the  adsor¬ 
bed  H2  molecules  are  in  the  range  of  0.88-0.95  A  in  the  three  cases. 
Analysis  of  the  charge  distributions  indicates  that  obvious  chemical 
bonds  are  also  formed  between  the  adsorbed  hydrogen  atoms. 
Therefore,  they  are  still  adsorbed  in  a  molecular  manner.  However, 
from  Table  1  we  can  see  that  the  average  H2  binding  energy  is 
always  larger  than  1.0  eV/H2.  This  is  too  large  to  be  reversible  for 
hydrogen  storage.  Thus,  V,  Cr  and  Mn  are  not  promising  adsorbates 
for  hydrogen  storage. 

The  adsorbed  Co,  Ni,  and  Cu  can  adsorb  two  H2  molecules  per 
metal,  and  their  configurations  are  shown  in  Fig.  2(b)— (d), 
respectively.  We  find  from  the  figures  that  Co  is  still  adsorbed  on 
the  H-site,  while  Ni  and  Cu  change  their  adsorption  sites  after  the 


adsorption  of  H2  molecules.  Ni  goes  to  the  S-site  and  Cu  goes  to  the 
D-site  without  any  energy  barrier.  Through  going  to  a  new 
adsorption  site,  the  charge  carried  by  the  adsorbate  atom  increases. 
For  example,  the  charge  carried  by  Cu  adsorbed  on  the  H-site  is 
only  +0.02e,  while  it  increases  to  +0.27e  after  going  to  the  D-site. 
The  increased  charge  polarizes  the  H2  molecules  and  decreases  the 
total  energy  of  the  system.  The  calculated  binding  energies  are 
1.26  eV/H2,  0.92  eV/H2  and  0.60  eV/H2  for  Co,  Ni  and  Cu  cases, 
respectively.  However,  the  gravimetric  density  is  only  2.8  wt%  in 
the  three  cases,  too  low  for  applications. 

3.2.  Hydrogen  storage  in  C6H6M2  complexes 

In  this  section,  we  investigate  the  hydrogen  storage  in  CeHeM2 
complexes,  where  two  metal  atoms  are  adsorbed  on  the  benzene 
surfaces,  one  on  each  side.  Binding  properties  of  various  metals  and 
their  hydrogen  storage  properties  are  listed  in  Table  2.  From  the 
table  we  can  see  metals  such  as  Be,  Na,  Mg,  K  and  Zn  still  cannot  be 
adsorbed  on  the  benzene  surfaces,  just  as  they  do  in  n  =  1  case.  Of 
course,  they  cannot  be  used  as  adsorbates  for  hydrogen  storage  in 
this  paper.  C6HeAl2  assumes  an  A-A  configuration,  where  two  Al 
atoms  occupy  two  A-sites  on  each  side  of  the  benzene  surfaces. 
However,  when  trying  to  add  H2  molecules  to  the  complex,  they 
always  escape.  Similar  things  happen  to  C6H6Ga2.  Therefore,  Al  and 
Ga  are  not  suitable  adsorbates  for  hydrogen  storage  too.  The 
adsorbed  Li  and  Ca  will  occupy  the  H-sites  to  form  H-H  configu¬ 
rations,  and  each  of  them  can  adsorb  four  H2  molecules.  Thus, 
CeH6Li2  and  CeH6Ca2  can  adsorb  up  to  eight  H2  molecules,  respec¬ 
tively.  We  note  that  each  adsorbed  Ca  atom  in  graphene  case  can 
also  adsorb  four  H2  molecules  [14].  However,  each  adsorbed  Li  atom 
can  adsorb  five  and  three  H2  molecules  in  C6o  [7]  and  MOF-5  [19] 
cases,  respectively.  Atomic  geometries  are  similar  for  CeH6Li2  and 
CeH6Ca2,  so  only  the  later  one  is  shown  as  an  example  in  Fig.  3(a). 
From  the  figure  we  can  see  that  H2  molecules  distribute  symmet¬ 
rically  around  the  adsorbate  atoms,  just  as  they  do  in  n  =  1  case.  The 
CeH6Li2  +  8H2  configuration  yields  an  extremely  high  H2  gravi¬ 
metric  density  of  14.8  wt%,  much  higher  than  requirement  (9  wt%) 
of  the  U.  S.  DOE  target  by  2015  [34].  More  importantly,  the  binding 
energy  is  0.29  eV/H2,  suitable  for  reversible  hydrogen  storage.  This 
is  larger  than  that  in  the  Li-decorated  C6o  [7]  and  MOF-5  [19]  cases, 
which  is  about  0.08  and  0.12  eV/H2,  respectively.  We  think  the 
different  Li  binding  energies  may  be  related  with  the  differences  of 
the  electronegativity  of  the  six-carbon  rings  in  the  three  materials. 
The  H2  density  in  CeH6Ca2  +  8H2  is  9.2  wt%  and  the  binding  energy 
is  0.45  eV/H2,  also  higher  than  the  2015  target  and  suitable  for 
reversible  hydrogen  storage.  Therefore,  Li  and  Ca  are  very  prom¬ 
ising  adsorbates  for  hydrogen  storage  when  benzene  molecules  are 
used  as  adsorbents. 

Two  transition  metal  atoms  can  also  adsorbed  on  the  benzene 
surfaces  to  form  double-decorated  complexes  and  then  adsorb  H2 


Fig.  1.  Atomic  geometries  of  (a)  C6H6Ca  +  4H2,  and  (b)  C6H6Ti  +  4H2  in  a  ball  and  stick  style;  (c)  charge  distributions  of  C6H6Ti  +  4H2  configuration  in  a  space-filling  style  with  an 
isosurface  value  of  0.08e/A3. 
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Fig.  2.  Atomic  geometries  of  (a)  C6H6Mn  +  3H2)  (b)  C6H6Co  +  2H2,  (c)  C6H6Ni  +  4H2)  and  (d)  C6H6Cu  +  2H2  configurations  in  a  bail  and  stick  style. 


molecules.  From  Table  2  we  can  see  the  bader  charge  on  the  metals 
decreases  from  Sc  to  Zn.  This  is  because  the  electronegativity 
increases  in  the  same  order.  Due  to  the  decrease  of  the  charge,  the 
metals  polarize  fewer  and  fewer  H2  molecules.  Sc  and  Ti  each  can 
adsorb  four  H2  molecules,  and  V  to  Fe  each  can  adsorb  three  ones, 
while  Co  to  Cu  each  can  adsorb  two  ones  only.  Sc  and  Ti  atoms  are 
likely  to  be  adsorbed  on  the  H-sites  to  form  C6FI6M2  complexes  in 
H— FI  or  sandwich-like  configurations.  After  the  adsorption  of  H2 
molecules,  they  still  have  similar  atomic  geometries,  and  we  show 
the  CeFIeTh  +  8H2  configuration  in  Fig.  3(b)  as  an  example.  The  FI2 
density  in  Sc  case  is  8.7  wt%,  very  close  to  the  U.  S.  DOE  target  by 
2015  [34],  and  the  binding  energy  (0.57  eV/ H2)  is  also  suitable  for 
reversible  hydrogen  storage.  Though  H2  density  in  Ti  case  is  also 
high,  the  binding  energy  is  somewhat  high  for  reversible  hydrogen 
storage.  In  V,  Cr,  Mn,  and  Fe  cases,  C6FI6M2  complexes  also  occur  in 
H-FI  configurations  and  each  complex  can  adsorb  six  H2  molecules 
at  most.  After  FI2  adsorption,  the  C6FI6M2  frameworks  are  almost 
unchanged.  Once  again,  due  to  their  similarities,  we  only  show  the 
atomic  geometry  of  CeFl6Cr2  +  6H2  in  Fig.  3(c)  as  an  example.  The 
adsorption  of  the  six  H2  molecules  results  in  H2  gravimetric 
densities  from  6.3  wt%  to  5.9  wt%.  Particularly,  the  density  in  Mn 
case  (see  Table  2)  just  equals  the  requirement  (6  wt%)  of  the  U.  S. 
DOE  target  by  2010  [34].  Flowever,  the  binding  energy  in  the  four 
cases  is  always  larger  than  1.1  eV/ H2,  too  large  to  be  reversible  for 
hydrogen  storage. 

In  Co,  Ni,  and  Cu  cases,  each  C6FI6M2  can  adsorb  four  H2  mole¬ 
cules,  and  the  corresponding  atomic  geometries  are  shown  in 
Fig.  3(d)-(f),  respectively.  The  three  metals  are  all  stable  at  the 
Fl-sites  before  H2  adsorption,  but  this  is  not  the  case  after  that.  Now 
we  examine  them  one  by  one.  From  Fig.  3(d)  we  can  see  after  H2 
adsorption  one  Co  atom  relaxes  to  the  upper  site  of  the  trapezoid 
(P4-site)  composed  of  the  four  neighboring  carbon  atoms,  while  the 


Table  2 

Metal  and  H2  binding  properties  in  n  =  2  (C6H6M2)  case. 


Metal 

M-site 

Eh(  M) 
(eV) 

M-charge 

Nmaxi  H2) 

Eb(  H2) 
(eV) 

H2  weight 
(wt%) 

Li 

H-H 

0.62 

+0.88e 

8 

0.29 

14.8 

Be 

- 

-0.20 

- 

- 

- 

- 

Na 

- 

-0.06 

- 

- 

- 

- 

Mg 

H-H 

0.09 

+0.05e 

- 

- 

- 

Al 

A-A 

1.62 

+0.34e 

- 

- 

- 

K 

- 

-0.09 

- 

- 

- 

- 

Ca 

H-H 

1.04 

+1.19e 

8 

0.45 

9.2 

Sc 

H-H 

3.11 

+0.93e 

8 

0.57 

8.7 

Ti 

H-H 

3.41 

+0.84e 

8 

0.90 

8.4 

V 

H-H 

2.48 

+0.69e 

6 

1.24 

6.3 

Cr 

H-H 

1.20 

+0.67e 

6 

1.20 

6.2 

Mn 

H-H 

1.51 

+0.47e 

6 

1.12 

6.0 

Fe 

H-H 

2.23 

+0.48e 

6 

1.18 

5.9 

Co 

P4-P3 

2.77 

+0.34e 

4 

1.62 

3.9 

Ni 

D-D 

2.83 

+0.1 6e 

4 

1.40 

3.9 

Cu 

S-D 

1.06 

+0.05e 

4 

0.68 

3.8 

Zn 

H-H 

0.16 

-O.Ole 

- 

- 

- 

Ga 

A-A 

1.34 

+0.42e 

- 

- 

- 

other  one  relaxes  to  the  upper  site  of  the  triangular  (P3-site) 
composed  of  the  three  neighboring  carbon  atoms.  Finally,  each 
adsorbed  Co  atom  adsorbs  two  H2  molecules,  yielding  a  H2  uptake 
of  3.9  wt%.  Flowever,  the  average  H2  binding  energy  is  rather  large 
(1.62  eV/ H2),  indicating  the  Co-benzene  complex  is  not  suitable  for 
hydrogen  storage.  In  Ni  case,  both  Ni  atoms  relax  to  the  D-sites  after 
FI2  adsorption,  one  on  each  side  of  the  benzene  surfaces,  and  each 
metal  adsorbs  two  FI2  molecules.  Again,  due  to  the  large  FI2  binding 
energy  and  the  low  gravimetric  density,  Ni-decorated  benzene 
complex  is  also  not  suitable  for  hydrogen  storage.  Each  Cu  can  also 
adsorb  two  FI2  molecules,  and  thus  the  C6H6CU2  complex  can 
adsorb  four  FI2  molecules  totally.  Flowever,  after  FI2  adsorption,  one 
Cu  goes  to  the  S-site  and  the  other  goes  to  the  D-site  as  shown  in 
Fig.  3(f).  The  binding  energy  is  0.68  eV/H2,  suitable  for  reversible 
hydrogen  storage,  but  the  gravimetric  density  is  3.8  wt%,  lower 
than  the  requirement  of  applications. 

3.3.  Dimerization  of  metals  on  the  benzene  surfaces 

From  above  discussions,  we  find  that  at  least  Li,  Ca  and  Sc  are 
promising  adsorbates  when  benzene  molecules  are  used  as 
adsorbents  for  hydrogen  storage.  Flowever,  there  are  still  questions. 
Will  these  metals  dimerize?  If  they  do  dimerize,  what  effects  will 
they  have  on  the  hydrogen  storage  capability?  To  answer  these 
questions,  we  will  consider  the  dimerization  of  various  metals  on 
the  benzene  surfaces  as  well  as  their  hydrogen  storage  properties  in 
this  section. 

Dimerization  studies  are  only  focused  on  those  metals  that  can 
be  strongly  adsorbed  on  the  benzene  surfaces  as  discussed  in  above 
sections.  We  first  place  two  atoms  of  one  metal  on  the  same  side  of 
the  benzene  surfaces  and  then  relax  the  structure.  After  full 
relaxations,  we  find  that  the  distance  between  the  two  metal  atoms 
is  always  small.  For  example,  the  equilibrium  Li — Li  distance  on  the 
benzene  surface  is  2.38  A,  much  smaller  than  the  calculated  bond 
length  of  2.76  A  of  an  isolated  Li2  dimer.  It  indicates  that  a  Li2  dimer 
is  formed  on  the  benzene  surface.  Similar  thing  happens  to  other 
metals,  so  they  will  form  their  corresponding  dimers  on  the 
benzene  surfaces  too.  Flowever,  to  one  metal  which  form  is  more 
stable?  Is  the  dimer  form  or  the  isolated  form?  We  calculate  and 
compare  the  total  energies  of  the  two  forms  to  every  metal,  finding 
that  all  the  metals  in  this  work  except  Ca  have  lower  energy  in  the 
dimer  form  than  in  the  isolated  form.  Therefore,  Ca  atoms  will 
prefer  to  be  isolated  while  other  metals  will  prefer  to  dimerize  on 
the  benzene  surfaces. 

We  first  focus  our  attention  on  the  metals  that  are  more  likely  to 
dimerize.  In  general,  dimerization  of  metals  is  less  favorable  for 
hydrogen  storage.  For  example,  dimerization  of  Li  atoms  on  C60 
surface  reduces  the  weight  percentage  of  hydrogen  storage  [11  ].  We 
have  known  from  above  discussions  that  Li  and  Sc  are  promising 
adsorbates,  so  we  only  consider  their  hydrogen  storage  capability 
after  dimerization.  Li2-  and  Sc2-dimer  benzene  complexes  are 
geometrically  similar,  so  we  only  shown  the  later  structure  in 
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Fig.  3.  Atomic  geometries  of  (a)  C6H6Ca2  +  8H2)  (b)  C6H6Ti2  +  8H2,  (c)  C6H6Cr2  +  6H2)  (d)  C6H6Co2  +  4H2,  (e)  C6H6Ni2  +  4H2,  and  (f)  C6H6Cu2  +  4H2  configurations  in  a  ball  and  stick 
style. 


Fig.  4(a).  From  the  figure  we  can  see  that  the  formation  of  the  metal 
dimer  slightly  distorts  the  benzene  geometry.  By  adding  H2  mole¬ 
cules  one  by  one  to  the  adsorbed  and  dimerized  metals,  we  find 
that  l\2  and  SC2  can  finally  adsorb  four  H2  molecules  per  dimer  at 
most.  The  geometry  after  FI2  adsorption  in  the  Sc2-dimer  case  is 
shown  in  Fig.  4(b)  as  an  example  again.  Flere  the  H2  densities  are 
8.0  wt%  and  4.6  wt%  for  U2-  and  Sc2-dimer  benzene  complexes, 
respectively,  remaining  high  or  moderate.  The  calculated  binding 
energies  are  0.25  eV/ H2  and  0.48  eV/ H2,  respectively,  still  suitable 
for  reversible  hydrogen  storage.  Therefore,  Li  and  Sc  are  promising 
adsorbates  for  hydrogen  storage  compared  with  other  metals, 
though  they  will  dimerize  on  the  benzene  surfaces. 

Finally,  we  would  like  to  do  further  discussions  on  Ca  metal.  As 
just  discussed,  two  Ca  atoms  will  prefer  to  be  isolated  on  the 
benzene  surfaces,  one  on  each  side,  forming  an  H — H  or  a  sand¬ 
wich-like  configuration  as  discussed  in  Section  3.2.  According  to 
our  calculations,  this  configuration  is  about  0.54  eV  lower  than  the 


Fig.  4.  Atomic  geometries  of  Sc2-dimer  benzene  complex  in  a  ball  and  stick  style  (a) 
before  and  (b)  after  H2  adsorption. 


dimer  one.  In  order  to  test  the  stability  of  the  CeFIeCan  complexes 
further,  first-principles  molecular  dynamics  simulations  are  per¬ 
formed  on  CeFIeCa  +  4H2  and  CeFl6Ca2  +  8FI2  systems.  The  simula¬ 
tions  are  performed  at  300  °C  with  a  time  step  of  lfs  for  lOps. 
Results  indicate  that  though  some  FI2  molecules  escape  from  the 
metals,  the  complex  frameworks  are  always  stable.  Combined  with 
the  high  H2  densities  and  the  suitable  H2  binding  energies,  we  can 
judge  that  Ca  is  the  most  promising  adsorbate  for  reversible 
hydrogen  storage. 

4.  Conclusions 

In  summary,  we  have  performed  first-principles  studies  on  the 
stability  and  hydrogen  storage  properties  of  various  metal-benzene 
complexes.  We  find  Ca  is  the  most  promising  metal  when  forming 
complexes  with  benzene  for  hydrogen  storage.  Besides,  Li-  and 
Sc-decorated  benzene  complexes  can  also  store  a  large  number  of 
H2  molecules  with  suitable  H2  binding  energies,  so  they  are  also 
promising  for  reversible  hydrogen  storage.  Other  metals  we  studied 
are  not  suitable  for  reversible  hydrogen  storage.  We  think  that  this 
work  will  be  interesting  for  the  studies  on  hydrogen  storage  in 
materials  containing  six-carbon  rings  such  as  C60,  CNTs,  MOFs  and 
so  on,  especially  for  the  selection  of  metal  adsorbates.  Finally  we 
hope  this  study  could  motivate  further  experimental  and  theoret¬ 
ical  studies  on  hydrogen  storage  in  metal-benzene  complexes  and 
benzene  ring  contained  materials. 
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